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We investigate numerically the coupled diffusion-advective type field equations originating from 
the canonical phase space approach to the noisy Burgers equation or the equivalent Kardar-Parisi- 
Zhang equation in one spatial dimension. The equations support stable right hand and left hand 
solitons and in the low viscosity limit a long-lived soliton pair excitation. We find that two identi- 
cal pair excitations scatter transparently subject to a size dependent phase shift and that identical 
solitons scatter on a static soliton transparently without a phase shift. The soliton pair excitation 
and the scattering configurations are interpreted in terms of growing step and nucleation events in 
the interface growth profile. In the asymmetrical case the soliton scattering modes are unstable 
presumably toward multi soliton production and extended diffusive modes, signalling the general 
non-integrability of the coupled field equations. Finally, we have shown that growing steps per- 
form anomalous random walk with dynamic exponent z — 3/2 and that the nucleation of a tip is 
stochastically suppressed with respect to plateau formation. 
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I. INTRODUCTION 

There is a continuing interest in the strong coupling 
aspects of stochastically driven nonequilibrium systems. 
The phenomena in question are ubiquitous and comprise 
turbulence in fluids, interface and growth problems, and 
chemical and biological systems. 

In this context the noisy Burgers equation or the equiv- 
alent Kardar-Parisi-Zhang (KPZ) equation, describing 
the nonequilibrium growth of a noise-driven interface, 
provide a simple continuum model of an open driven non- 
linear system exhibiting scaling and pattern formation. 

In one dimension, which is our concern here, the noisy 
Burgers equation for the local slope, u(x,t) = ^7h(x,t), 
of a growing interface has the form 



du o . ^ , 
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The height profile (in a comoving frame) h(x, t) is then 
governed by the equivalent KPZ equation 



(3) 



In (Q) and (||) v is the damping or viscosity characteriz- 
ing the linear diffusive term, A a coupling strength for the 
nonlinear mode coupling or growth term, and n a Gaus- 
sian white noise, driving the system into a statistically 
stationary state. The noise is correlated according to 
(||) and characterized by the strength A. Moreover, the 
Burgers equation is invariant under the slope-dependent 
Galilean transformation 



i.e., the interface is superimposed with a constant slope 
in a moving frame. 

The Burgers equation (|l|) and its KPZ equivalent in 
one and higher dimensions have been the subject of in- 
tense scrutiny in recent years owing to its paradigmatic 
significance within the field theory of nonequilibrium sys- 
tems §-{l5). 

In a series of papers the one dimensional case defined 
by (Q) and (^) has been analyzed in an attempt to un- 
cover the physical mechanisms underlying the pattern 
formation and scaling behavior. Emphasizing that the 
noise strength A constitutes the relevant nonperturba- 
tive parameter that is driving the system into a statisti- 
cally stationary state, the method was initially based on 
a weak noise saddle point approximation to the Martin- 
Siggia-Rose functional formulation of the noisy Burgers 
equation fl6|-|i9|]. This work was a continuation of ear- 
lier work based on the mapping of a solid-on-solid model 
onto a continuum spin model |p0|j . More recently the 
functional approach has been superseded by a canonical 
phase space method deriving from the canonical structure 
of the Fokker-Planck equation associated with the Burg- 
ers equation Below we briefly summarize these 
findings. 

The functional or the equivalent phase space approach 
valid in the weak noise limit, A — + 0, replaces the stochas- 
tic Langevin-type Burgers equation (|l|) by coupled deter- 
ministic diffusion-advection type mean field equations, 
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for the slope u(x, t) and a canonically conjugate noise 
field p(x,t), replacing the stochastic noise r\. The field 
equations bear the same relation to the Fokker-Planck 
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equations as the classical equations of motion bear to 
the Schrodinger equation in the semi-classical WKB ap- 
proximation [ p6| . 

To justify the weak noise limit we recall the anal- 
ogy with the WKB approximation in quantum mechan- 
ics which, owing to its nonperturbative character, cap- 
tures features like bound states and tunneling ampli- 
tudes, which are generally inaccessible to perturbation 
theory. Therefore, we anticipate that the present weak 
noise approach to the Burgers equation also accounts cor- 
rectly, at least in a qualitative sense, for the stochas- 
tic properties even at larger noise strength. However, 
there may be an upper threshold value beyond which 
the system may enter a new stochastic or kinetic phase. 
In the one dimensional case discussed here the scaling 
behavior is controlled by a single strong coupling fixed 
point which can be accessed by the present weak noise 
approach. In two and higher dimension a dynamic renor- 
malization group analysis predicts a kinetic phase tran- 
sition at a critical noise strength (or coupling strength) 
and the weak noise approach presumably fails. 

The equations (^|) and (^) derive from a principle of 
least action characterized by an action S(u' — » u",t) as- 
sociated with an orbit u'(x) — > u"{x) traversed in time t 



/•*»« / du 
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J0,u' 

with Hamiltonian density 



dt 



H=p[ vX7 2 u + XuVu - -\7 2 p 
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The action is of central importance in the present ap- 
proach and serves as a weight function for the noise- 
driven nonequilibrium configurations in much the same 
manner as the energy E in the Boltzmann factor 
exp(— (3E) for equilibrium systems, where (3 is the in- 
verse temperature. The dynamical action in fact replaces 
the energy in the context of the dynamics of stochastic 
nonequilibrium systems governed by a generic Langevin 
equation driven by Gaussian white noise. The action pro- 
vides a methodological approach and yields access to the 
time dependent and stationary probability distributions, 



P(v! — > u" , t) oc exp 



S(u' -> u", t) 



P s t(u") = lim P{u — > u",t) , 

t — >oc 
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and associated moments, e.g., the stationary slope corre- 
lations 



u(xt)u(00)) = 

u"(x)u'(0)P(u' 



u",t)P at (u') 
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The canonical formulation associates the conserved en- 
ergy E (following from time translation invariance), the 



E = J dxH , 
II = J dx uVp , 
M = dxu . 



conserved momentum II (from space translation invari- 
ance), and the conserved area M (from the Burgers equa- 
tion with conserved noise): 



(12) 
(13) 
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The field equations (^) and (0) determine orbits in a 
canonical up phase space where the dynamical issue in 
determining S and thus P is to find an orbit from u' 
to u" in time t, p being a slaved variable. Note that 
unlike dynamical system theory we are not considering 
the asymptotic properties of a given orbit. In general 
the orbits in phase space lie on the manifolds determined 
by the constants of motion E, II and M. Here the zero 
energy manifold E = plays a special role in defining 
the stationary state. For vanishing or periodic boundary 
conditions for the slope field the zero energy manifold is 
composed of the transient submanifolds p = and the 
stationary submanifold p = 2vu. The zero energy or- 
bits on the p — manifold correspond to solutions of 
the damped noiseless Burgers equation; the orbits on the 
p = 2vu are solutions to the undamped noiseless Burgers 
equation with negative damping, i.e., v replaced by —v. 
In the solvable linear case of the noise driven diffusion 
equation for A = 0, i.e., the Edwards- Wilkinson equation 
p8[ , a finite energy orbit from u' — > u" in time t migrates 
to the zero energy manifold in the limit t — > oo, yielding 
according to Eqs. (0) and ( |l0|) the stationary distribu- 
tion P st oc exp(— (^/A) J dxir). This distribution also 
holds in the Burgers case and is a generic result indepen- 
dent of A p8[| . Finally, in the long time limit an orbit 
from u' — > u" is attracted to the hyperbolic saddle point 
at the origin in phase space implying ergodic behavior 
in the stationary state. In Fig. 1 we have schematically 
depicted possible orbits in phase space. 

The field equations (||) and (^) admit nonlinear soliton 
or smoothed shock wave solutions which are, in the static 
case, of kink-like form, 



u\{x) = u tanh 



AH 
2v 



(15) 



Propagating solitons are subsequently generated by the 
Galilean boost (|J). Denoting the right and left bound- 
ary values by u+ and u_, respectively, the propagation 
velocity is given by 



u + + it_ = —2vj\ 



(16) 
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p=2uu, E=0 
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the numerical method designed in order to treat the in- 
herent instability. In Section 4 we present our numerical 
results for the scattering of two single solitons on a static 
soliton and the scattering of two soliton pairs. In Sec- 
tion 5 we discuss growth and nucleation associated with 
the modes investigated numerically. Section 6 is devoted 
to a summary of our results and a conclusion. 



FIG. 1. Generic behavior of the orbits in up phase space. 
Heavy lines indicate the zero energy manifold. The stationary 
saddle point (sp) is at the origin. The finite time orbit from 
u' to u" is attracted to the saddle point for t — > oo. 

The amplitude of the static soliton is u and the soli- 
ton is located at the origin. The right hand soliton for 
u > 0, i.e., the soliton with the larger right hand side 
boundary value, moves on the noiseless manifold p = 
and is also a solution of the damped (stable) noiseless 
Burgers equation for rj = 0. The noise-induced left hand 
soliton for u < 0, i.e., the soliton with the larger left hand 
side boundary value, is associated with the noisy mani- 
fold p — 2isu, and is a solution of the undamped (unsta- 
ble) noiseless Burgers equation with v replaced by — v. 
In addition the field equations also admit linear mode 
solutions superimposed as ripple modes on the solitons 
p0| . The ripple modes are superpositions of both decay- 
ing and growing components reflecting the noiseless and 
noisy manifolds p — and p — 2vu, respectively. The 
soliton mode induces a propagating component with ve- 
locity Xu in such a way that the right hand soliton acts 
like a sink and the left hand soliton as a source of linear 
modes. In the Edwards- Wilkinson limit for A — > the 
ripple modes become the usual diffusive modes (growing 
and decaying) of the driven stationary diffusion equation. 

The heuristic physical picture that emerges from our 
analysis is that of a many body formulation of the pattern 
formation of a growing interface in terms of a dilute gas 
of propagating solitons matched according to the soliton 
condition ( |l6| ) with superimposed linear ripple modes. In 
Fig. 2 we have depicted the basic soliton modes constitut- 
ing the building blocks in the representation of a growing 
interface. For further illustration we have shown in Fig. 3 
the slope field u, the corresponding height field h, and the 
noise field p for a 4-soliton configuration. 

In the present paper we embark on a numerical anal- 
ysis of the coupled field equations (||) and (g) with the 
purpose of investigating them in more detail and provide 
a numerical underpinning of the heuristic quasi-particle 
picture advanced in the work referred to above. The pa- 
per is organized in the following manner. In Section 2 
we discuss the soliton modes. In Section 3 we introduce 



FIG. 2. Slope field u and height profile h for the right 
hand and left hand moving kink solitons forming the 'quarks' 
in the description of a growing interface. 
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FIG. 3. 4-soliton representation of the slope field u, the 
noise field p, and the height field h. 
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II. SOLITON MODES 



The exact right and left hand soliton solutions of 
the field equations do not satisfy periodic or vanishing 
boundary conditions in the slope field u; the nonvanish- 
ing boundary values u + and U- in fact correspond to a 
deterministic current dissipated or generated at the soli- 
ton centers yielding permanent profile solutions |3lJ]. The 
kink solitons constitute the elementary building blocks or 
'quarks' in the present approach and the interface profile 
is then built up by matching solitons according to the 
matching condition (|l6|). 

The simplest mode satisfying periodic boundary con- 
ditions is the two-soliton or pair soliton configuration 



ui{x,t) = U\{x — vt — Xi) — i*i [x — vt — X2) 



(17) 



obtained by matching a right hand and a left hand soliton 
boosted to the velocity v = —\u. The two-soliton mode 
has amplitude 2u and size \xz —Xi\. The associated noise 
field vanishes for the right hand component and equals 
2vu for the left hand component; we thus have 

ps(x, t) = — 2vu\(x — vt — X2) for u>0, (18) 
Pi{x, t) = +2vu\{x — vt — x\) for u<0. (19) 

By inspection it is seen that the pair mode ( |l7j ) is an 
approximate solution to the field equations (||) and (||). 
The correction terms are of the type uVu and uVp re- 
ferring to the distinct components of 1*2 and "pi and thus 
correspond to local perturbations from a region of size 
v/\\u\ which is small in the low viscosity limit v — > 0. 
We assume that the correction can be treated within a 
linear stability analysis and thus gives rise to a linear 
mode propagating between the right hand and left hand 
solitons ]30|] . 

The pair mode thus forms a long-lived excitation or 
quasi-particle in the many body description of a grow- 
ing interface. Subject to periodic boundary conditions 
this mode corresponds to a simple growth situation. The 
propagation of the pair mode corresponds to the propa- 
gation of a step in the height field h. At each revolution 
of the pair mode the interface grows by a uniform layer of 
thickness 2u\x2 — x\\. In Fig. 4 we have depicted the pair 
mode in u, the associated noise field p, and the height 
profile h. 

Generally a growing interface, ignoring the superim- 
posed linear ripple modes, can at a given time instant 
be represented by a gas of matched left hand and right 
hand solitons as depicted in Fig. 3 in the 4-soliton case. 
A gas of pair solitons thus constitute a particular growth 
mode where the height profile between moving steps has 
horizontal segments. 
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FIG. 4. Slope field 112, the associated noise field pi, and 
the resulting height profile hi at time t = for a two-soliton 
configuration. 

The soliton picture also allows us easily to understand 
in what sense the right hand soliton acts like a drain and 
the left hand soliton as a source with respect to pertur- 
bations. Considering two pair solitons superimposed on 
the right and left horizontal parts of the static solitons 
@ it follows from © that for a right hand soliton per- 
turbations move toward the soliton center and for a left 
hand soliton perturbations move away from the soliton 
center. This mechanism also follows from the linear ana- 
lysis of ripple modes in |p0| . The mechanism is depicted 
in Fig. 5 



FIG. 5. Source and drain mechanism for the right hand 
and left hand solitons. The perturbation attracted and re- 
pelled by the soliton centers are modeled by pair solitons. 

Since we do not posses explicit propagating multi- 
soliton solutions of the field equations (||) and (^|) the 
problem of soliton collisions remains unresolved from an 
analytical point of view. Therefore we now turn to a 
numerical analysis of the problem. 
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III. NUMERICAL METHOD 

The coupled field equations (||) and (^|) are of the 
diffusion-advection type with the exception that the evo- 
lution of p is governed by a negative diffusion coefficient. 
Standard numerical methods designed to step the equa- 
tions forward in time fail because small perturbations 
with wavenumber k grow in time like exp(z^fc 2 i), so per- 
turbations with the largest possible k grow fastest, ren- 
dering hence the integration unstable. In order to circum- 
vent this problem we have designed a method to solve the 
equations iteratively starting with a trial solution for u 
and p in the space interval \x\ < L and time interval 
< t < T. At each iteration step we proceed in two 
sweeps. In the forward sweep we step only the equation 
for u forward in time using p(x, t) from the previous iter- 
ation step. In the backward sweep we step the equation 
for p backward in time, using u(x,t) from the forward 
sweep. In this manner the unwanted perturbations in p 
decrease exponentially as one moves backward in time. 

A drawback of this method is that we can specify ini- 
tial values only for u, and that we must instead specify p 
at the final time t = T. In the present context where we 
want to consider collisions of solitons this is not a serious 
problem, because here it is possible to guess the final so- 
lution. The numerical solution serves therefore mainly as 
a tool to check that a certain guess is actually a solution. 
Furthermore, this method allows us to calculate the pre- 
cise functional form of u and p during the collision, even 
if the initial guess around the time of the collision was 
actually wrong. 

We have solved the equations on a mesh with N x = 
1001 mesh points and N t = 6001 time steps. For L = 1/2 
we have a mesh spacing of Ax = 0.001. For both sweeps 
we use sixth order finite differences to calculate first 
and second derivatives and a third order Runge-Kutta 
scheme for the time integration (see, e.g., the appendix 
of Ref . |Q for these schemes) . 

In order to adequately resolve u and p at all times 
we must choose a suitable value of v. We found em- 
pirically that v = 0.005 gave good results, which is the 
value adopted in the following. For smaller values of v 
the u and p functions become only marginally resolved 
whereas for larger values of v the length of the time step 
is mostly controlled by the value of v rather than just 
the propagation speeds of the solitons. Empirically we 
found that the maximum time step that can be used is 
At = 5 x 10~ 5 for v = 0.005. In all cases we have chosen 
A = 1. 



IV. SOLITONS, PAIR SOLITONS, AND SOLITON 
COLLISIONS 



In choosing soliton configurations to be verified by the 
numerical method we have found that it is essential to 



satisfy the three conservation laws governing the dynam- 
ics of solitons, namely the conservation of energy jl^), 
the conservation of momentum ([l3]), and the conserva- 
tion of area (HJ) . 



A. Solitons and pair soliton 

We have numerically verified that the right and left 
hand solitons ( |l5| ) for u > and u < are solutions. 



By construction the two-soliton configuration (17) carries 
energy E 2 — (— 16/3)A^|u| 3 , momentum H 2 = — 4^u|u|, 
and for small v the area M 2 oc 2v\x\ —x 2 \ ■ We have shown 
that, for small v, a well-separated two-soliton mode, i.e. 
with \x\ — x 2 \ 3> v/Xu, is a long-lived excitation, hence 
justifying the heuristic argument. 

In order to lend support to the heuristic quasi-particlc 
picture based on the elementary kink solitons (quarks) 
and the composite pair soliton as the basic quasi-particle 
it is essential to consider soliton collisions. We have here 
considered the two symmetric cases: i) the collision of two 
propagating solitons with a static soliton and ii) the colli- 
sion of two pair solitons. In both cases the configurations 
are symmetric and the conservation laws are satisfied at 
all times including the collision regime. 



B. Three soliton collisions 

In the first case two propagating solitons moving in op- 
posite directions collide with a static soliton located at 
the center. The trial solution has the form 

u(x, t) = — sign(t) [u± [x + vt) + ui(x — vt) — 2«i(2x)] 

(20) 

p(x,t) = -2v[ui(x + vt) +ut(x-vt)] for t<0 (21) 
p(x,t) = -Avu-l{2x) for t>0 (22) 

and the height field 
h(x, t) = — sign(t) 

lop 



2v 
T 



cosh 


^(x-vt) 


cosh 


'^(x + vt) 


cosh 


X\u\ 
2u ^ 





(23) 



with velocity v = Xu. 

In this mode two left hand solitons with amplitude 2u 
propagate with equal and opposite velocities toward a 
static right hand soliton with amplitude 4m located at the 
center. During the collision the left hand solitons are ab- 
sorbed, the static right hand soliton flips over to a static 
left hand soliton and two right hand solitons emerge prop- 
agating away from the center with equal and opposite 
velocity. The solitons thus collide transparently with the 
static soliton, i.e., there is no reflection, and there is no 
phase shift associated with the collision. In terms of the 
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t = 0.00 



associated height profile this scattering situation corre- 
sponds to filling in a dip with subsequent nucleation of a 
growing tip. 

Energy and momentum are associated with the noise- 
induced left hand solitons moving on the noisy manifold 
p = 2vu. By inspection of ( f2p| ) it follows that the total 
energy E = — (32/3)^Au 3 , the total momentum II = 0, 
and the total area M = are conserved during the colli- 
sion. 

Choosing the amplitude u = 2 we have in Fig. 6 de- 
picted the numerical verification of the slope field u as a 
function of x for different values of t. In Figs. 7 and 8 
we have shown the associated noise field p and the height 
profile h as a function of x for the same values of t. In 
Fig. 9 we have shown a gray-scale representation of u in 
the xt plane. We notice that there is no phase shift as- 
sociated with the scattering process. Finally, in Fig. 10 
we have shown the convergence of energy, momentum, 
and area during the forward and backward time sweeps 
integrating numerically the field equations. 
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FIG. 7. The noise field p as a function of x for the same 
values of t as in Fig. 6. 
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FIG. 6. Three soliton collision: The slope field u as a 
function of x for different values of t. 



FIG. 8. The height profile h = J u Ax as a function of x 
for the same values of t as in Fig. 6. 
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C. Pair soliton collisions 



FIG. 9. Gray-scale representation of u in the xt plane 
showing the three soliton collision. Note the absence of a 
phase shift dugii^fhe collision. 
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FIG. 10. Evolution of E, II, and M during the first few 
iteration steps. The dashed line indicates the first iteration 
step and the solid line the last one. Dotted lines are for in- 
termediate steps. 



In the second case we consider the collision of two pair 
solitons of equal size and amplitude. The trial solution 
propagating with velocity v = Xu has the form 

u(x, t) = — [u\(x — vt + x\) ~ ui(x — vt + X2)] 
+ [u\{x + vt — x 2 ) — ui(x + vt — Xi)] 



p(x, t) = —2v[ui{x — vt + Xi) + Ui{x + vt — Xi)] 
for < vt < X2, 
u(x, t) = — [ui(x — vt + X\) — ui(x)] 

+ [u 1 (x) - Ul(x + Vt- .Tl)] 

p(x, t) — —2v[ui(x — vt + Xi) + Ui(x + vt — Xi)] 
for X2 < vt < xi, 

u(x, t) — +[ui(x + vt — Xi) — ui(x)] 
— [u\{x) — Ui(x — vt + Xi)] 
p(x, t) = — 2v[ui(x) + ui(x)] 

for x\ < vt < 2x\ — X2, and 

u(x,t) = +[u\(x + vt — x\) — u\(x + vt — 2x\ + X2)] 
— [u\(x — vt + 2x\ — X2) — u\(x — vt + x\)} 



(24) 
(25) 



(26) 
(27) 



(28) 
(29) 



p(x, t) = —2v\u,\(x + vt + 2x\ — X2) 
+ui(a; — vt + 2x\ — X2)] 



(30) 



(31) 



for 2x\ — X2 < vt. 

In this mode two pair solitons of amplitude 2u prop- 
agate with equal and opposite velocities toward one an- 
other. The two leading kink solitons merge to a static 
soliton and the two trailing kinks are absorbed. Subse- 
quently, the static right hand soliton flips over to a static 
left hand soliton and the two pair solitons re-emerge. 
Analyzing the collision it follows that the scattering of 
pair solitons is transparent and accompanied by a phase 
shift in space equal to the soliton size \x2 — x\ \ or, equiv- 
alently, a time delay \x2 — x\\jv. In terms of the asso- 
ciated height profile the scattering situation corresponds 
to filling in a trough due to two colliding steps and the 
subsequent nucleation of a growing plateau. 

By inspection it again follows that the total energy 
E = — (32/3)^Au 3 , the total momentum II = 0, and the 
total area M = are conserved during collision. 

Choosing the amplitude 2u and the kink positions 
x\ = 0.25 and X2 — 0.15 we have shown in Fig. 11 the 
numerical verification of the slope field a as a function 
of x for different values of t. In Figs. 12 and 13 we have 
shown the associated noise field p and the height profile 
h as a function of x for the same values of t. In Fig. 14 
we have shown a gray-scale representation of u in the xt 
plane. We notice the phase shift engendered during the 
transparent collision. Finally, in Fig. 15 we have shown 
the convergence of energy, momentum, and area during 
the forward and backward time sweeps in the numerical 
integration. 
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FIG. 11. Pair soliton collision: The slope field it as a 
function of x for different values of t. 
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FIG. 12. The noise field p as a function of x for the same 
values of t as in Fig. 11. 



FIG. 13. The height profile h = J u Ax as a function of x 
for the same values of t as in Fig. 11. 




FIG. 14. Gray-scale representation of u in the xt plane 
showing the pair soliton collision. Note the occurrence of a 
phase shift during the collision. 
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V. GROWTH AND NUCLEATION 

Since we have achieved numerical justification of three 
specific dynamical soliton configurations, namely i) the 
pair soliton, ii) the collision of two solitons with a static 
soliton, and iii) the collision of two pair solitons, we can 
proceed to draw some simple conclusions based on the 
general framework discussed in Sect. 1. There are two 
levels of description: the stochastic Langevin level and 
the deterministic Fokker-Planck or equation of motion 
level. On the Fokker-Planck level yielding the canoni- 
cal field equations (||) and (||) the growth of the inter- 
face is interpreted in terms of a gas of propagating soli- 
tons (and diffusive modes). The stochastic description 
on the Langevin level is then established in the weak 
noise limit A — » by computing the action S associated 
with a particular dynamical mode and subsequently de- 
duce the probability distribution according to (^), i.e., 
P oc exp(— S/A). This procedure is completely equiva- 
lent to the WKB limit of quantum mechanics. Here the 
wavefunction $ and thus the probabilistic interpretation 
is given by $ oc exp(iS/h), where S is the action asso- 
ciated with the classical motion ^(|. Note that unlike 
quantum mechanics there is no phase interference in the 
stochastic nonequilibrium case. 



A. The pair soliton: anomalous diffusion 

The long-lived pair soliton ( \L7\ ) has size £ = \x\ — x%\, 
amplitude 2u, and propagates with velocity v = — Alt. 
During a revolution in a system of size L with peri- 
odic boundary conditions the height field increases with 
a layer of thickness 2u£. Since the system is traversed 
in time t = L/v the integrated growth velocity is given 
by 2Xu 2 £/L which for a single pair of fixed size vanishes 
in the thermodynamic limit. On the other hand, the lo- 
cal growth velocity dh/dt is given by 2Xu 2 = (A/2)(V/i) 2 
which is consistent with the averaged KPZ equation (0) 
in the stationary state. 

The stochastic properties of the pair soliton growth 
mode is also easily elucidated by noting that the ac- 
tion associated with the pair mode is given by S — 
(4/3)i/A|u| 3 i. Denoting the center of mass of the pair 
mode by x = (x\ + x-i)j2 we have u = vj A = x/tX and 
we obtain using (^|) the transition probability 

P(x,t)(xexp(-^-£^^j , (32) 

for the 'random walk' of independent pair solitons or 
steps in the height profile. Comparing ( |32"|) with the 
distribution for 'ordinary' random walk originating from 
the Langevin equation dx/dt = i] , (r/rj) (t) = AS(t), 
P(x,t) oc exp(— x 2 /2At), we observe that the growth 
mode performs anomalous diffusion. The distribution 
( |32| ) also implies the soliton mean square displacement, 
assuming pairs of the same average size, 

(a 2) (t)oc ^_j (33) 

with dynamic exponent z = 3/2, identical to the dynamic 
exponent defining the KPZ universality class. This result 
should be contrasted with the mean square displacement 
(x 2 ) oc At 2 ' z , z = 2, for ordinary random walk. The 
growth modes thus perform superdiffusion. 

B. Soliton collisions: nucleation 

The three soliton scattering case shown in Figs. 6, 7, 
and 8 corresponds to the filling in of a dip and the sub- 
sequent nucleation of a tip. The incoming solitons have 
amplitude 2u and move with velocities v = ±Xu. Denot- 
ing the distance of a soliton from the center by d, the 
height change from the bottom of the dip to the top of 
the tip is Ah = Aud. The duration of the collision process 
is At = 2d/v and we obtain for the local growth velocity 
Ah/ At = (X/2)(2u) 2 = (A/2)(V/i) 2 , consistent with the 
KPZ equation. 

As regards the stochastic properties the action prior 
to the collision is associated with the incoming left hand 
solitons and given by (8/3)vXu 3 t. After the collision the 
static left hand soliton of amplitude 4u carries the action 







(32/3)v\uH. There is thus an increase of action in con- 
nection with the tip formation yielding, according to @ , 
a reduced probability, i.e., a statistical suppression of tip 
nucleation. 

The case of pair soliton scattering shown in Figs. 11, 
12, and 13 correspond to the filling in of a trough and the 
subsequent formation of a growing plateau. The incom- 
ing pair solitons have amplitude 2u and propagate with 
velocities v — ±Xu. Denoting the pair size by d the height 
change from the trough to the plateau is Ah = Aud, the 
duration of the collision At = 2d/v, and we obtain again 
Ah/ At = (A/2)(V7i), consistent with the KPZ equation. 

Before and after the collision the action is now un- 
changed and equal to (8/3)v\u 3 t. As in the three soliton 
case this action comes from the two left hand solitons. 
However, during the collision in a time interval At = d/v 
the action jumps to the value (32/3)vXu 3 t as the static 
left hand soliton of amplitude Au is temporarily formed 
before the re-emergence of the pair solitons. Comparing 
this situation with the tip formation in the three soliton 
scattering case we conclude that tip formation is statisti- 
cally suppressed with respect to the formation of a grow- 
ing plateau. Generally the formation of tips is 'expensive' 
since a large change in slope yields a large change in ac- 
tion. 



VI. SUMMARY AND CONCLUSION 

In the present paper we have numerically investigated 
the coupled diffusion-advection type field equation origi- 
nating from the canonical phase space approach applied 



to the noisy Burgers equation or the equivalent KPZ 
equation in one spatial dimension. We have shown that 
the pair soliton mode in the slope field corresponding to 
moving steps in the height field forms a long-lived ex- 
citation. We have furthermore investigated two special 
scattering scenarios, namely the collision of two identical 
moving solitons with a static soliton and the collision of 
two identical pair solitons. They correspond in the height 
field respectively to the nucleation of a growing tip and 
to the formation of a growing plateau. Finally, we have 
applied the canonical phase space approach in order to es- 
timate the stochastic aspects of the above configurations 
and found i) that a step in the height field performs a ran- 
dom walk with dynamical exponent z = 3/2 correspond- 
ing to averaged superdiffusion and ii) that tip formation 
is stochastically suppressed with respect to plateau for- 
mation. 

As discussed above the inherently unstable structure 
of the field equation makes a direct integration forward 
in time inaccessible and we can thus not establish solu- 
tions as an initial value problem and discuss the equations 
generically. Consequently, we are limited to numerically 
check trial solutions representing a variety of scattering 
situations. So far we have only been able to verify the 
symmetric cases of three soliton and soliton pair colli- 
sions. In order to extend the present numerical approach 
and thus provide substance to the heuristic quasi-particle 
representation of a growing interface it is clearly of inter- 
est to design more involved trial solutions. Alternatively, 
a completely different approach to generate solutions is 
called for. 
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